Vibration Mode Induced Shapiro Steps and Backaction in Josephson Junctions 
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A model of a superconducting tunnel junction coupled to a mechanical oscillator is studied at zero 
temperature in the case of linear coupling between the oscillator and tunneling electrons. It is found 
that the Josephson current flowing between two superconductors is modulated by the motion of the 
oscillator. Coupling to harmonic oscillator produces additional Shapiro steps in I — V characteristic 
of Josephson junction whose position is tuned by the frequency of the vibration mode. We also find 
a new velocity dependent term originating from the backaction of the ac Josephson current. This 
term is periodic in time and vanishes at zero bias voltage. 
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The coupling of the charge carriers to vibrational 
modes and localized spins in electronic devices has been 
a subject of intense investigation recently. Vibrational 
modes and spins possess dynamic internal degrees of free- 
dom, much unlike static impurities or defects. As a con- 
sequence, they affect the electronic dynamics in these 
devices. Interesting I-V characteristics (i.e., peaks in 
differential tunneling conductance) in molecular electron- 
ics 0,0, 0,0, HI El may indicate strong influence from the 
electronic-vibrational coupling. A step structure (rather 
than peak structure) in the differential tunneling con- 
ductance has also been observed in the STM-based in- 
clastic tunneling spectroscopy around a local vibrational 
mode on surfaces |jj . The vibrational effects on the con- 
ductance of molecular quantum dots were also examined 

nam in iia. 

Spin detection and manipulation is crucial in spintron- 
ics and quantum information processing. The modu- 
lation of the current by precessing spins may be used 
to detect and manipulate single spins in electron-spin- 
resonance-scanning tunneling microscopy (ESR-STM) 
technique [Tl ITi l ll H Il7j . Theoretically, it is known [13, 
LLii Uii LLlI that spin-orbit interactions and the polariza- 
tion of injected electrons are important in modulating the 
tunneling current. Theses studies focused on the tunnel- 
ing between two normal metals, where only the single 
particle process is involved in transport. Recently, we 
studied a precessing spin coupled through a local direct 
exchange interaction to the tunneling electrons between 
two superconductors [5^ . 

We now consider, for the first time, the Josephson ef- 
fect in a superconducting tunneling junction coupled to a 
mechanical oscillator in the tunnel barrier. Both effects 
of the oscillator motion on the tunnel current and the 
tunneling electrons on the oscillator dynamics are stud- 
ied. To our knowledge, none of these effects have been ad- 
dressed before. We find that: (i) In the tunnel junction, 
the Josephson current is modulated by the motion of the 
oscillator; the Fourier spectrum of the Josephson current 
exhibits peaks at frequency luj ±Wo, 2wj, and ujj ± 2loq, 
where ujj = 2eV and luq are the Josephson frequency 
and the vibrational mode of the oscillator respectively. 




V 



FIG. 1: An oscillator in the tunneling barrier coupled to two 
superconducting leads. The motion of the oscillator mod- 
ulates the electron tunneling between the superconductors. 
The chemical potentials on both sides differ by the applied 
voltage, fiL — fiR = eV. 



These additional peaks are the result of coupling of ac 
Josephson current at wj and oscillation at usq. (ii) Elec- 
tron tunneling through the junction leads to a novel time 
dependent change in oscillator energy. If no voltage bias 
is applied across the junction (i.e., dc Josephson effect in 
equilibrium), the oscillator energy is time independent. 
When a nonzero voltage bias applied, the oscillator will 
display time-dependent energy variations. 

Our model is illustrated in Fig. ^ It consists of two 
ideal superconducting leads coupled to each other by a 
tunneling barrier, which contains a mechanical oscillator. 
The Hamiltonian 



H = Hr 



Hr 



H 7 



(1) 



The first two terms are, respectively, the Hamiltonians 
for electrons in the left and right superconducting leads 
of the tunnel junction, 



H L(R) - Ek(p);<7 e k(p) c k(p),<T C k(p),a 



-£k(p)[ A L(fl)C k ( P )T ct - 



k(-p)i 



H.c] , 



(2) 



where wc denote the electron creation (annihilation) op- 
erators in the left, right (L,R) leads by (cker) and 



"per 



(c pcr ) respectively k (p) are momenta and a is the spin 
index, ewp),<7 an d ^l(r) are , respectively, the single par- 
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ticle energies of conduction electrons, and the pair po- 
tential (the gap function) in the leads. Without loss of 
generality, we assume that the superconductors are of 
a spin-singlet s-wave pairing symmetry, and consider the 
Josephson tunneling at zero temperature. The third term 
in Eq. Q depicts the tunneling between the supercon- 
ductors: 



Hi 



kp c kCT C pcr 



H.C. 



k,p;er 



where the tunneling matrix elements T kp transfer elec- 
trons through an insulating barrier. When a local vibra- 
tional mode is embedded into the tunneling barrier then, 
in the linear coupling regime, 



(4) 



where a describes the coupling between the tunneling 
electrons and vibrational mode. The quantity u is the 
displacement operator for the oscillator. 

Henceforth, we denote by M and K the mass and 
spring constant of the mechanical oscillator. As the 
energy associated with the vibrational mode, ujq = 
y/K/M ~ 10 _1 - 10~ 6 eV is much smaller than the 
typical electronic energy on the order of 1 eV, the me- 
chanical oscillation is very slow as compared to the time 
scale of electronic processes. This allows us to apply the 
Born-Oppenheimer approximation to treat the electronic 
degrees of freedom as if the local oscillator is static at 
every instantaneous location. We will treat the dynam- 
ics of the mechanical oscillator including the back action 
from the tunneling electrons. 

Whenever a voltage bias is applied across the junction, 
the Josephson current 



Ij(t) 



dt>[e ieV ( t+t '){[A(t),A(t')]_) 

OO 

e- ieV( - t+t '\[A\t),A\t')]-)\ , (5) 



(that is, a = 0), which is found to be, 



4\eV) = e£ 
k, P 



| A |2| T (0) |2 



1 



eV + E k + E p 



eV-E k 
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(3) with |A| is the superconducting energy gap and -E"k(p) 



\/( e k(p) - E F ) 2 + |A| 2 . The amplitude 



T S (eV) 
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k P 



(E u 
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r(0) 



An order of magnitude estimate gives Tsujj/Jg 
{eV/\A\) 2 « 1. 

From Eq. ©, we may construct the coupling modu- 
lated part of the Josephson junction energy 



Hj = Ej{l + au) 2 [l-cos(ujjt)} 

+— a{l + cm)— sm(u>jt) 
2e at 



(9) 



where Ej = /2e. The derivative of Hj with respect 
to the phase yields the supercurrent in Eq.©. Eq.@ 
captures the Josephson backaction effect- it vividly il- 
lustrates how the electronic degrees of freedom influence 
the mechanical oscillator of Fig.l coupled to them. We 
may obtain this result by also employing other standard 
techniques. For instance, integrating out the fermionic 
degrees of freedom, we obtain an effective action S e f / for 
u(t) in Matsubara time. When the standard mechanical 
part of the oscillator is added, the action 

S = \ fdr jdr' \T^\ 2 (l + au(T))(l + au(r'))K(r,r') 



1 f d 

where the operator A(t) - £ k , p:CT T kp dl a (t)c pa (t). +- I dr[Ku 2 + M(^) 2 } 

Here c k(p)(T (t) = e iif M«)* Ck(p)(7 e-^W* with K L(R) = 

H L (R) - VL(R)N L(R) and N L(R) = £k(p),<r c{, p)(T c k(p)(T . The non-locality in time present 



(10) 



The unequal chemical potentials of the two superconduc- 
tors lead to a voltage bias /xl — /i_r? = eV . Hereafter, we 
set H = 1. A little algebra yields: 



Ij(t) = J ( s Q \eV)[l + au\ 2 sm(u;jt) 

. du 

+T s (eV)(l + au)a— cos(ujjt) 



(6) 



The non-locality in time present in the first term reflects 
how electronic correlations between different times influ- 
ence the oscillator coordinates and lead to the Josephson 
backaction effect. Here, K{t,t') = Fl(t — t')F r {t' — 
t) cxp[eV r (r + r')], with F the Matsubara-Gorkov func- 
tion. As the resulting classical action is not time trans- 
lationally invariant, energy is not conserved. Within the 
Born-Oppenheimer approximation (u(t') ~ u(r) + (t' — 
T)du/dr), the equation of motion 



Here we assume that the two superconductors are iden- 
tical and set the constant phase difference between them 
4>q = 0. The Josephson frequency is given by loj = 2eV. 
The quantity Jg is the amplitude of the Josephson cur- 
rent in the absence of coupling to the vibrational mode 



= 



ss 

Su(t 



M—+Ku 
dt z 



\Tg!\ 2 J dr'K(r, r')[l + au(r) + a^(r' - r)], (11) 
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leading to 



, r d 2 u , . du rr _ . . 
M w+lf s(t)^ + Ku = F{t) 



(12) 



Here, the driving force 



F(t) = -2aEj(l + au)y.-[l + ^^]cos(ujjt)\ , (13) 

with T s ujj/AeEj ~ (eV/|A|) 2 < 1, and the time- 
dependent energy non-conserving 



7 S (t) = -(a 2 r s /e)sm(wjt) 



(14) 



These relations are also very transparent within the real 
time approach of Eq. (jjj, where the total oscillator 
Hamiltonian H osc = Hj + P 2 /2M + Ku 2 /2. A more 1 
detailed Kcldysh contour analysis yields the same results 
|23| . In Ea. p3fl . the first, linear in u term, may, alterna- 
tively, be lumped into the spring moduli K and regarded 
as a Josephson stiffness- a shift of the spring constant 
resulting from electronic correlations. The oscillatory 
part of the driving force and the time dependence of the 
Josepshon backaction generated 7s(£) (Ea.^l4jl) lead to 
interesting experimental consequences which we will now 
elaborate on. Equations JfjJ), (|12fl . (|13(l and (|14|l consti- 
tute the central result of this work. 

Josephson Bakaction. As shown by Eq. I)14|) and 
Eq. JSJl, the velocity coefficient (75 (t)) originates from 
the coupling of the mechanical oscillator to the tunneling 
electrons. It is quadratically proportional to the coupling 
constant a, which is similar to the case of normal metal 
tunnel junctions Q. This term has two novel features: 
(i) 75 depends on the voltage bias. At zero voltage bias 
(the dc setting), 7s vanishes since T$ is zero. 7s is finite 
only when a finite voltage bias is applied across the junc- 
tion (i.e., ac case). In the low voltage limit (eV -C |A|), 
75 is linearly proportional to the voltage bias, (ii) Once 
a finite voltage bias is applied, 7s is also a periodic func- 
tion of time with the Josephson frequency uij. Both of 
these features are absent in the normal metal tunnel junc- 
tions (24) . These properties are unique to the coupling 
of the mechanical oscillator to the superconductors. In 
the normal metal, there is no quasiparticle energy gap on 
the Fermi surface. Matters become far richer when the 
oscillator is coupled to the superconductor. On the one 
hand, in the superconductor, there exists an energy gap 
on the Fermi surface and the quasiparticles are depleted 
below this energy. This leads to the quenching of the 
single particle tunneling channel; no contribution to the 
dissipation of the oscillator due to normal quasiparticles 
is possible. On the other hand, due to macroscopic quan- 
tum coherence in the superconductor, Cooper pairs can 
tunnel through the barrier between the two superconduc- 
tors with a probability comparable to that of single par- 
ticle tunneling in a normal metal junction. When a static 
voltage bias is applied, the tunneling of a single Cooper 
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FIG. 2: Absolute magnitude of the Fourier transform Ij(uj) of 
the Josephson current as given by Eq. J I Oil for the vibration 
mode u>o = OAljj (a) and 1.4 (b). The other parameters 
values are taken to be: ce = 0.1, and K = 0.6. 



pair requires an energy 2eV to overcome the potential 
barrier. The energy carried by Cooper pair upon tunnel- 
ing can be transferred to oscillator. For zero bias voltage, 
the tunneling pairs do not acquire/lose energy (the effec- 
tive electronic only action is time independent, no effec- 
tive external sources are present, and energy is preserved 
at all times). In the ac setting, smtujt is odd under time 
reversal and, as a consequence, the term sinwji^ is al- 
lowed in the effective Hamiltonian Eq. 10. 

Shapiro steps. To calculate the Josephson current, we 
need to solve Eq. 1|12|) for the displacement field u(t). 
In the weak coupling limit and in view of the fact that 
Ts is much smaller than J^, the main physics can be 
captured by neglecting the damping terms and the aTs 
and a 2 terms in the driving force. In this limit, 



u(t) = uocos(o;ot) + 



2aE, 



M{u>l - w>) 
and the Josephson current 



cos(wji) 



2aEj 



K 



(15) 



Ij(t) = ^ (1 — -75-) sm(ujjt) + 2<5cos(woi) sm(uj jt) 



a 2 cos 2 (coot) sm(ujt) > ,(16) 



4a 2 
K 

2a 2 sin(2wji) 



where a = auo and K = Kuq/Ej. Equation (|16|) demon- 
strates clearly that the Josephson current not only oscil- 
lates with time with a frequency uij, but is also modu- 
lated by the vibrational mode of the mechanical oscillator 
with a frequency loq- In Fig. El wc plot the absolute mag- 
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nitude of the Fourier transform of the Josephson current 
given by Eq. i|16fl for various values of the vibration mode 
frequency ojq. The spectrum shows a main peak at the 
frequency ujj. In addition, the coupling of the mechani- 
cal oscillator and tunneling electrons generates new side 
peaks at frequencies ojj ± loq, ojj ± 2ujq, and 2u>j. The 
intensity of these peaks is proportional to the coupling 
constant. Note that the peak at 2uij originates from the 
second term in u(t) given by Eq. fTK|l . which is a direct 
manifestation of the feedback effect from the Josephson 
tunneling. Our calculation implies that a dc component 
arises if the voltage bias has one of the Shapiro step values 
uq and 2wo. When higher-order effects are taken into ac- 
count, the equation-of- motion for the oscillator can only 
be solved numerically. The main conclusions presented 
here remain qualitatively unchanged. 

Conclusion. We studied the Josephson junction cou- 
pled to a mechanical oscillator between its two supercon- 
ducting leads. We found that the Josephson current flow- 
ing between two spin-singlet pairing superconductors is 
modulated by the motion of the oscillator. The coupling 
of an oscillator of cigcnfrequcncy luq to an ac junction 
of characteristic frequency ujj leads to beats. We find 
novel Shapiro steps induced at wj ± uj , ujj ± 2ujq, and 
2lu,j. This differs from the case of a precessing spin j22j . 
where due to the sum rule of tunneling through different 
spin channels, the Josephson current is not modulated if 
the spin-orbit coupling mechanism does not exist. The 
electrons tunneling through the superconducting junc- 
tion lead to a novel non-energy conserving effect. If a 
voltage bias is applied, this time dependent effect arises 
from the back-action of the supercurrent on the oscilla- 
tor dynamics. As far as we know, no measurements of 
Josephson current through a vibrational mode between 
two superconductors have been reported yet. Recent 
progress in molecular electronics and nanomechanical 
resonators jl^,|2^| holds great promise in attaching single 
molecules to superconducting leads, or even tune the tun- 
nel barrier of the superconducting junctions by a mechan- 
ical cantilever. Our predictions are, potentially, within 
the realm of present technology. Another possible exper- 
iment concerns atomically sharp superconducting tip in 
low t emp erature STM in both the quasiparticle tunnel- 
ing |27l | and Josephson tunneling regimes [2^| (coined 
"Josephson STM" or JSTM ^) on conventional super- 
conductors. It is very interesting to extend the JSTM 
technology by using a superconducting tip to study the 
Josephson current in the vicinity of a local vibrational 
mode on the superconducting surface, which may pro- 
vide a new detection technique. 
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